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Abstract This chapter reviews the nature of turbulence in the Galactic interstellar 
medium (ISM) and its connections to the star formation (SF) process. The ISM is 
turbulent, magnetized, self-gravitating, and is subject to heating and cooling pro- 
cesses that control its thermodynamic behavior, causing it to behave approximately 
isobarically, in spite of spanning several orders of magnitude in density and temper- 
ature. The turbulence in the warm and hot ionized components of the ISM appears 
to be trans- or subsonic, and thus to behave nearly incompressibly. However, the 
neutral warm and cold components are highly compressible, as a consequence of 
both thermal instability in the atomic gas and of moderately-to-strongly supersonic 
motions in the roughly isothermal cold atomic and molecular components. Within 
this context, we discuss: i) the production and statistical distribution of turbulent 
density fluctuations in both isothermal and polytropic media; ii) the nature of the 
clumps produced by thermal instability, noting that, contrary to classical ideas, they 
in general accrete mass from their environment in spite of exhibiting sharp dis- 
continuities at their boundaries; iii) the density-magnetic field correlation (and, at 
low densities, lack thereof) in turbulent density fluctuations, as a consequence of 
the superposition of the different wave modes in the turbulent flow; iv) the evolu- 
tion of the mass-to-magnetic flux ratio (MFR) in density fluctuations as they are 
built up by dynamic compressions; v) the formation of cold, dense clouds aided by 
thermal instability, in both the hydrodynamic (HD) and the magnetohydrodynamic 
(MHD) cases; vi) the expectation that star-forming molecular clouds are likely to be 
undergoing global gravitational contraction, rather than being near equilibrium, as 
generally believed, and vii) the regulation of the star formation rate (SFR) in such 
gravitationally contracting clouds by stellar feedback which, rather than keeping the 
clouds from collapsing, evaporates and diperses them as while they collapse. 
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1 Introduction 



The interstellar medium (ISM) of our galaxy (the Milky Way, or simply, the Galaxy) 
is mixture of gas, dust, cosmic rays, and magnetic fields that occupy the volume in- 
between stars. The gasesous component, with a total mass ~ 1O 1O M0, may be in 
either ionized, neutral atomic or neutral molecular forms, spanning a huge range 
of densities and temperatures, from the so-called hot ionized medium (HIM), with 
densities « ~ 10~ 2 cm~ 3 and temperatures T ~ 10 6 K, through the warm ionized 
and neutral (atomic) media (WIM and WNM, respectively), with n ~ 0.3 cm~ 3 and 
T ~ 10 4 K and the cold neutral (atomic) medium (CNM, n ~ 30 cm -3 , T ~ 100 
K), to the giant molecular clouds (GMCs, n > 100 cm" 3 and T ~ 10-20 K). GMCs 
can span several tens of parsecs across, and, in turn, contain plenty of substructure, 
which is commonly classified into clouds (n ~ 10 3 cm -3 , size scales L of a few 
parsecs), clumps (n ~ 10 4 cm -3 , L ~ 1 pc), and cores (n > 10 s cm" 3 , 0.1 pc). 
It is worth noting that the temper ature of most molecular gas is remarkably uniform, 



10-30 K (e.g.. lFerrierell2001h . 



Moreover, the ISM is most certainly turbulent, as typical estimates of the Reynolds 
number (j? e ) within it are v ery large. For example, in the cold ISM, R e ~ 10 5 -10 7 
dElmegreen & Scaloll2004 §4.1). This is mostly due to the very large spatial scales 
involved in interstellar flows. Because the ISM's temperature varies so much from 
one type of region to another, so does the sound speed, and therefore the turbu- 
lent velocity fluctuations are often moderately or even strongly supersonic (e.g., 
Heiles & Troland , 2003 ; Elmegreen & Scaloi 2004 , and references therein). In these 



cases, the flow is significantly compressible, inducing large-amplitude (nonlinear) 
density fluctuations. The de nsity enhancements thus formed constitute den se clouds 



and their substructure (e.g., Sasao . 1973 : Ballesteros-Paredes et al. . 1999ah 



In addition to being turbulent, the ISM is subject to a number of additional phys- 
ical processes, such as gravitational forces exerted by the stellar and dark matter 
components as well as by its own self-gravity, magnetic fields, cooling by radia- 
tive microscopic processes, and radiative heating due both to nearby stellar sources 
as well as to diffuse background radiative fields. It is within this complex and dy- 
namical medium that stars are formed by the gravitational collapse of certain gas 
parcels. 

In this chapter, we focus on the interaction between turbulence, the effects of 
radiative heating and cooling, which effectively enhance the compressibility of the 
flow, the gas' self-gravity, and magnetic fields. Their complex interactions have a 
direct effect on the star formation process. The plan of the chapter is as follows: in 
§S]we briefly recall the effects that the net heating and cooling have on the effective 
equation of state of the flow and, in the case of thermally unstable flows, on its ten- 
dency to spontaneously segregate in distinct phases. Next, in S|3]we discuss a few 
basic notions about turbulence and the turbulent production of density fluctuations 
in both the hydrodynamic (HD) and magnetohydrodynamic cases, to then discuss, 
in §2] the evolution and properties of clouds and clumps formed by turbulence in 
multiphase media. In $5] we discuss the likely nature of turbulence in the diffuse 
(warm and hot) parts of the ISM, as well as in the dense, cold atomic and molecular 
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clouds, suggesting that in the latter, at least those in the process of forming stars, 
the velocity field may be dominated by a gravitational contraction. Next, in Sj6]we 
discuss the regulation of star-formation (SF) in gravitationally contracting molec- 
ular clouds (MCs), in particular whether it is accomplished by magnetic support, 
turbulence, or stellar feedback, and how. Finally, in S|7]we conclude with a summary 
and some final remarks. 



2 ISM Thermodynamics: Thermal Instability 



The ISM extends essentially over the entire disk of the Galaxy and, when consid- 
ering a certain dense subregion of it, such as a cloud or cloud complex, it is nec- 
essary to realize that any such subregion constitutes an open system, whose inter- 
actions with its environment need to be taken into account. A fundamental form 
of interaction with the surroundings, besides dynamical interactions, is through the 
exchange of heat, which occurs mostly due to heating by the U V background radia- 
tion produced by distant massive stars, local heating when nearby stellar sources 
(OB star ionization heating and supernova explosions), cosmic ray heating, and 
cooling by thermal and line emission from du st and gas, respectively (see, e.g., 



Dalg arno & McCravLll972tlWolfire et al.Ml995l) 



Globally, and as a first approximation, the ISM is roughly isobaric, as illustrated 
in Fig. 03 As can be seen there, most types of regions, either dilute or dense, lie 
within an order of magnitude from a thermal pressure P ~ 3000 K cm -3 . The largest 
deviations from this pressure uniformity are found in HII regions, which are the 
ionized regions around massive stars due to the star's UV radiation, and molecular 
clouds, which, as we shall see in ^5.41 are probably pressurized by gravitational 
compression. 

The peculiar thermodynamic behavior of the ISM is due to the functional forms 
of the radiative heating and cooling functions acting on it, which depend on the 
density, temperature, and chemical composition of the gas. T he left panel of Fig. [2] 
show s the temperature dependence of the cooling function A ( Dalgarno & McCravT 
19721) . One well-known crucial consequence of this general form of the cooling is 
that the atomic medium is thermally unstable dFieldl 119651) in the density range 



1 n £ 10 cm (corresponding to 5000 £ T ^ 30 K), meaning that the medium 
tends to spontaneously segregate into two stable phases, one warm and diffuse, with 
n ~ 0.3 cm~ 3 and T ~ 8000 K, and the other cold and dense, with n ~ 30 cm~ 3 and 

2500 K cm" 3 dField et al " 



80 K, both at a pressure P/k < 



1969; Wolfire et al 



1995, see also the reviews by Meerson 1996; Vazquez-Semadeni et al. 2003; 
Vazquez-Semadeni 2012), as illustrated in the right panel of Fig. [2] The cold gas 
is expected to form small clumps, since the fastest growing mode of the instabil- 
ity occurs at vanishingly small scales in the absence of thermal conductivit y, or at 
scales ~ 0.1 pc for the estimated thermal conductivity of the ISM (see, e.g. jField , 
1965tlAudit & Hennebellel 120051) . Because the atomic gas in the ISM has two stable 
phases, it is often referred to as a thermally bistable medium. 
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Fig. 1 Thermal pressure in various types of interstellar regions. The points labeled coronal corre- 
spond essentially to what we refer to as the HIM in the text; intercloud regions refer to th e WIM 
and W NM; diffuse, to CNM clouds, and dark, globule and molecular to molecular gas. From lMversI 
Clll). 



It is important to note that, even if the medium is not thermally unstable, the 
balance between heating and cooling implies a certain functional dependence of 
Peq(p), which is often approximated by a polytropi c law of the form P e q K p" (e.g., 
Elmegreeni 1199 It IVazquez-Semadeni et al. 1 19961) . where y e is the effective poly- 
tropic exponent. In general, y e is not the ratio of specific heats for the gas in this 
case, but rather a free parameter that depends on the functional forms of A and F. 
The isobaric mode of TI corresponds to y e < 0. A flow is sometimes said to be softer 
as the parameter y e becomes smaller. 



3 Compressible Polytropic MHD Turbulence 

3.1 Equations 

In the previous section we have discussed thermal aspects of the ISM, whose main 
dynamical effect is the segregation of the medium into the cold and warm phases. 
Let us now discuss dynamics. As was mentioned in SjT] the ISM is in general tur- 
bulent and magnetized, and therefore it is necessary to understand the interplay be- 
tween turbulence, magnetic fields, an d the effects of the net cool i ng (nA — F), which 
affects the compressibility of the gas ( Vazquez-Semadeni et al. , 19961) . 
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Fig. 2 Left: Thermal-equilibrium pressure P eq as a function of number density for "standard" 
conditions of metallicity and background UV radi ation for the atomic medium. The horizontal 
axis gives log 1() (/t/cm 3 ). From lWolfire et all jl995j) . Right: Temperature dependence of the cool- 
ing function. The labels ind icate values of the ionization fraction (per number) of the gas. From 
iDalgarno & McCravl fl972h . 



The dynamics of the ISM are governed by the fluid equations, complemented by 
self-gravity, the heating and coolin g terms in th e energy equation, and the equation 
of magnetic flux conservation (e.g.. lShulll992l) : 



dp 



du 
~dt 
de 
~dt 



= -pV-M, 

VP 



— + u ■ Vm = V<p + v 

P 



V(V-m) 



1 



(1) 



(Vx5)x B, (2) 



dB 



■ u-Ve= -{y- \)eV -u+r-nA 
|Vx(Bxu) = —V x (77 V x B) + V x 
V 2 cp = 47iGp, 



B 



4Jry c p n Pi 



4-np 

(3) 

, (4) 
(5) 



where p = ji m mnn is the mass density, [l m is the mean particle mass, mn is the 
hydrogen mass, u is the velocity vector, e is the internal energy per unit mass, B 
is the magnetic field strength, <p is the gravitational potential, v is the kinematic 
viscosity, 7] is the electrical resistivity, and y c is the collisional coupling constant 
between neutrals and ions in a partially ionized medium. Equation ([]]) represents 
mass conservation, and is also known as the continuity equation. Equation (O is the 
momentum conservation, or Navier-Stokes equation per unit mass, with additional 
source terms representing the gravitational force V<p/p and the Lorentz force. In 
turn, the gravitational potential is given by Poisson 's equation, eq. (0. Equation © 
represents the conservation of internal energy, with F being the heating function 
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and A the cooling function. The combination nA — F is usually referred to as the 
net cooling. Finally, eq. represents the conservation of magnetic flux (see below). 
Equations (d])-© are to be solved simultaneously, given some initial and boundary 
conditions. 

A brief discussion of the various terms in the momentum and flux conservation 
equations is in order. In eq. (0, the second term on the left is known as the advective 
term, and represents the transport of /-momentum by the j component of the veloc- 
ity, where i and j represent any two components of the velocity. It is responsible 
for mixing. The pressure gradient term (first term on the right-hand side [RHS]) in 
general acts to counteract pressure, and therefore density, gradients across the flow. 
The term in the brackets on the RHS, the viscous term, being of a diffusive nature, 
tends to erase velocity gradients, thus tending to produce a uniform flow. Finally, 
the last term on the RHS is the Lorentz force. 

On the other hand, eq. (01, assuming 77 = (i.e., zero electrical resistivity, 
or equivalently, infinite conductivity) and y c — > °° (i.e., perfect coupling between 
neutrals and ions), implies that the magnetic flux <P through a Lagrangian cross- 
sectional area A, given by 



remains constant in time as the area moves with the flow. This is the property known 
as flux freezing, and implies that the gas can slide freely along field lines, but drags 
the field lines with it when it moves perpendicularly to them. Note that this condition 
is often over-interpreted as to imply that the magnetic and the density fields must 
be correlated, but this is an erroneous notion. Only motions perpendicular to the 
magnetic field lines produce a correlation between the two fields, while motions 
parallel to the lines leave the magnetic field unaffected, while the density field can 
fluctuate freely. We discuss this at more length in 8 13.41 

The first term on the RHS of eq.|4]represents dissipation of the magnetic flux by 
electrical resistivity, and gives rise to the phenomenon of reconnection of field lines 
(see, e.g., the book by Shu 1992, and the review by Lazarian 2012). The second term 
on the RHS of eq. (j4]i represents ambipolar diffusion (AD), the deviation from the 
perfect flux-freezing condition that occurs for the neutral particles in the flow due to 
their slippage with respect to the ions in a partially ionized medium. We will further 
discuss the role of AD in the process of star formation (SF) in ^3.4.3l and l5.4.2l . 



3.2 Governing N on-Dimensional Parameters 




(6) 



Turbulence develops in a flow when the ratio of the advective term to the viscous 
term becomes very large. That is, 
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where R e is the Reynolds number, U and L are characteristic velocity and length 
scales for the flow, and G denotes "order of magnitude". This condition implies that 
the mixing action of the advective term overwhelms the velocity-smoothing action 
of the viscous term. 

On the other hand, noting that the advective and pressure gradient terms con- 
tribute comparably to the production of density fluctuations, we can write 



0{u-Vu) 

Wp/p) 



Ui r 
L 



AP 
LP 



Ap 
P 



U 2 



Mi 



= mI 



Ap_ 
P 



(8) 
(9) 



where M S = U /c s is the sonic Mach number, and we have made the approximation 
that AP I Ap <~~> c\, where c s is the sound speed. Equation (0 then implies that strong 
compressibility requires M s ^> 1. Conversely, flows with M s < 1 behave incom- 
pressibly, even if they are gaseous. Such is the case, for example, of the Earth's at- 
mosphere. In the incompressible limit, p = est., and thus eq. (]]) reduces to V ■ u = 0. 
Note, however, that the requirement M s ^> 1 for strong compressibility applies for 
flows that behave nearly isothermally, for which the approximation AP j Ap ~ c\ is 
valid, while "softer" ( cf. jj2]> flows have much larger d ensity jumps at a given Mach 
number. For example. IVazquez-Semadeni et al. (1996) showed that polytropic flows 
of the form P<x with y e — s- 0, have density jumps of the order e Ms . 

A trivial, but often overlooked, fact is that, in order to produce a density enhance- 
ment in a certain region of the flow, the velocity at that point must have a negative 
divergence (i.e., a convergence), as can be seen by rewriting eq. ([T]) as 



dp 
dt 



— = — pV- u, 



(10) 



where d/dt = d /dt + u- Vis the total, material, or Lagrangian derivative. However, 
it is quite common to encounter in the literature discussions of pre-existing density 
enhancements ("clumps") in hydrostatic equilibrium. But it should be kept in mind 
that these can only exist in multi-phase media, where a dilute, warm phase can have 
the same pressure as a denser, but colder, clump. But even in this case, time, formation 
of that clump must have initially involved the convergence of the flow towards the 
cloud, and the hydrostatic situation is applicable in the limit of very long times after 
the formation of the clump, when the convergence of the flow has subsided. 

Finally, two other important parameters determining the properties of a mag- 
netized flow are the Alfvenic Mach number, Ma = I// v a an d me plasma beta, 

= Pth/Pmag where va = B/y/Anp is the Alfven speed. Similarly to the non- 
magnetic case, large values of the Alfvenic Mach number are required in order to 
produce significant density fluctuations through compressions perpendicular to the 
magnetic field. However, it is important to note, as mentioned in £ 13-11 that compres- 
sions along the magnetic field lines are not opposed at all by magnetic forces. Note 
that, in the isothermal case, j3 = 2c\ jv\. 
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3.3 Production of Density Fluctuations. The Non-Magnetic Case 



As mentioned in the previous sections, strongly supersonic motions, or the abil- 
ity to cool rapidly, allow the production of large-amplitude density fluctuations in 
the flow. Note, however, that the nature of turbulent density fluctuations in a single- 
phase mediurrQ (such as, for example, a regular isothermal or adiabatic flow) is very 
different from that of the cloudlets formed by TI (cf. Sj2]). In a single-phase turbulent 
medium, turbulent density fluctuations must be transient, as a higher density gener- 
ally conveys a higher pressure^ and therefore the fluctuations must re-expand after 
the compression that produced them has subsided. 

Note that the above result includes the case with self-gravity, since in single- 
phase media, although hydrostatic equilibrium solutions do exist, they are gener- 
ally unstable. S pecifically, the singular isothermal sphere is known to be unsta- 
ble dShul 1977b. and non-singular co nfigurations such as the Bonnor-Ebert (BE) 
spheres ( EbertL 1955 : Bonnoi , 19561) need to be truncated so that the central-to- 
peripheral density ratio is smaller than a critical value ~ 14 in order to be stable. 
Such stable configurations, however, need to be confined by some means to prevent 
their expansion. But if such a warm phase is not available, the only way to confine 
the BE sphere is to continuously extend it to infinity, in which case the central- 
to-peripheral density ratio also te nds to infinity, and the configuration is unstable 



dVazquez-Semadeni et all 12005a) 



Instead, in multiphase flows, abrupt density variations may exist between dif- 
ferent phases even though they may be at roughly the same thermal pressure, and 
therefore, the dense clumps do not tend to re-expand. In this section we discuss 
the probability distribution of the density fluctuations, the nature of the resulting 
clumps and their interfaces with their environment, the correlation between the mag- 
netic and density fields, and the evolution of the mass-to-magnetic flux ratio as the 
clumps are assembled by turbulent compressions. 



3.3.1 The Probability Distribution of Density Fluctuations. The Non-Magnetic 

Case 

For astrophysical purposes it is important to determine the distribution of the den- 
sity fluctuations, as they may constitute, or at least provide the seeds for, what we 
normally refer to as "clouds" in the ISM. Because of the transient nature of turbu- 
lent density fluctuations in single-phase media, however, this distribution refers to a 



1 Thermodynamically, a phase is a re gion of space throughou t which all physical properties of a 
material are essentially uniform (e.g.. lModell & Reid LI 1974b . A phase transition is a boundary 
that separates physically distinct phases, which diffe r in most thermodynamic variables except one 
(often the pressure). See I Vazquez-Semadenil j2009h for a discussion on the nature of phases and 
phase transitions in the ISM. 

2 An exception would b e a so-ca l led B urgers' flow, which is characterized by the absence of the 
pressure gradient term iBurgersl 1 19741) . and can be thought of as the transitional regime from 
thermal stability to instability. 
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time-stationary population of fluctuations, although the fluctuations themselves will 
appear and disappear on timescales that are short compared to the time over which 
the distribution is considered. 

The probability density distribution (PDF) of the density field in turbulent isother- 



mal flows was initially investigated through numerical simulations. IVazquez-Semadeni 



(1994) found that, in the isothermal case, the PDF posesses a logno rmal form. A the- 

ory for the emergence of this functional form was later proposed by lPassot & Vazquez-Semadeni 
(119981) . in which the production of density fluctuations was assumed to arise from a 
succession of compressive or expansive waves, each one acting on the value of the 
density left by the previous one. Because the medium contains a unique distribution 
of (compressible) velocity fluctuations, and because the density jumps in isother- 
mal flow depend only on Mach number but not on the local density, the density 
fluctuations belong all to a unique distribution as well, yet each one can be consid- 
ered independent of the others if the global time scales c onsidered are much lo nger 
than the autocorrelation time of the velocity divergence ( Blaisdell et al. , ll993l) . Fi- 
nally, because the density jumps are multiplicative in the density (cf. eq. |9), then 
they are additive in s = In p. Under these conditions, the Central Limit Theorem can 
be invoked for the increments in s, implying that s will be normally distributed. In 

consequence, p will have a lognormal PDF. 

In addition, iPassot & Vazquez-Semadeni (1998) also argued that the variance of 
the density fluctuations should scale linearly with Ms, a suggestion that has b een 
investigated further by other groups dPadoan et al l ll997tlFederrath et allEoOsh . In 
particular, using numerical simulations of compressible turbulence driven by either 
solenoidal (or "vortical") or compressible (or "potential") forces, the latter authors 
proposed that the variance of s is given by 



Oi = ln(l +£M S 2 ), 



(11) 



where b is a constant whose value depends on the nature of the forcing, taking the 
extreme values of b = 1/3 for purely solenoidal forcing, and b = 1 for purely com- 
pressible forc ing. The lognormal density PDF for the one-dimensional, isothermal 
simulations o flPassot & Vazquez- Semadenildl998l). with its dependence on M s , i s il- 
lustrated in the left panel of Fig. [3] Finally, Passot & Vazquez-Semadenil(ll998i see 
also Padoan & Nordlund 1999) also investigated the case where the flow behaves as 
a polytrope (cf. i|2]i with arbitrary values of y e , by noting that in this case the sound 
speed is not constant, but rather depends on the density as c s p^ -1 )/ 2 , implying 
that the local Mach number of a fluid parcel now depends on the local density be- 
sides its dependence on the value of the flow velocity . Introducing this dependence 
of M s on p in the expression for the lognormal PDF, Passot & Vazquez-Semadenil 
(119981) concluded that the density PDF should develop a power-law tail, at high 
densities when % < 1, and at low densities when y e > 1. This result was then con- 
firmed by numerical simulations of polytropic turbulent flows (Fig. [3] right panel). 
Physically, the cause for the deviation of the PDF from the lognormal shape is that, 
for 7e > 1, the sound speed increases with increasing density, and therefore, high- 
density regions can re-expand and disappear quickly, while "voids", with a lower 
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Fig. 3 Left: Lognormal density PDFs for isothermal one-dimensional simulations at various Mach 
numbers, indicated by the labels. The independent variable is s = In p. Right: Density PDFs for 
polytropic case s (i.e., with P °e p 7e ), with effective p olytropic exponent y e = 0.3 (top) and y e = 1.7 
(bottom). From lPassot & Vazq uez-Sema denil jl998h . 



sound speed, last for long times. For y e < 1, the sound speed decreases for increas- 
ing density, and the behavior is reversed: large-amplitude density peaks have lower 
sound speeds and therefore last for longer times, while the voids have higher sound 
speeds and disappear quickly. The resulting topology of the density field is illus- 
trated in the one-dimensional case in Fig. [4] 



3.3.2 The Nature of Turbuent Clumps 

The Ambiguity of Clump Boundaries and Masses 

The clumps produced as turbulent density fluctuations are precisely that: fluctua- 
tions in a continuum. Besides, the re is in general a mass and energ y flux through 



any boundary we choose to define ( IBallesteros-Paredes et al.L 1 1999al) . This is espe 



cially true in isothermal flows, where no transition from a diffuse phase to a dense 
one with the same pressure can occur, so that the only density discontinuities pos- 
sible are those produced by shocks. In this case, the density fluctuations produced 
by turbulent compressions (a transient event of elevated ram pressure) must always 
be transients, and must eventually re-expand or collapse (see ^3 .3b . In thermally 
bistable media, the "boundaries" are somewhat better defined, although there is still 
mass flux through them (see £14.2. 11 1. 

The elusiveness of the notion of clump boundaries implies that we must rethink 
some of our dearest notions about them. First of all, the mass of a cloud or clump 
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Fig. 4 Plots of s = lnp of the density field in one-dimensional simu l ations of polytropic flows 
with a resolution of 6144 grid points by iPassot & Vazquez-Semadenil 11998T) . Top: A simulation 
with 7 e = 0.5, exhibiting high-amplitude, narrow density peaks, due to the low sound speed at high 
density. Bottom: A simulation with y e = 1-5, exhiting low-amplitude, extended density peaks and 
deep "voids" due to the high values of the sound speed at high densities, and low values at low 
densities. 



is not well defined, and additionally must evolve in time. It is ill-defined because 
the clump boundary itself is. Several procedures exist for "extracting" clumps from 
observational maps or from numerical simulations. One of the most wid ely em- 



ployed algorithms for locating clumps is Clumpfind (W illiams et al.L 119941) . which 



works by locating local peaks in the field being examined and then following the 
field down its gradient until another clump profile is met, at which point an arbitrary 
boundary is defined between the two. However, this procedure, by construction, is 
uncapable of recognizing "hierarchical", or "nested", structures, where one coherent 
"parent" clump contains other equally coherent "daughter" ones. Moreover, not sur- 
prisingly, it has been shown that the clump sets obtained from application of this al- 
gorithm depend sens itively on the parameters chosen for the definition of the clumps 



dPineda et al.L |2009). Conversely, a technique that, by definition, is capable of de- 



tecting "parent" structures, is based on "structure trees", or "dendrograms" (e.g., 
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Houlahan & Scald 1992 ; Rosolowskv et al. , 2008 ). which works by thresholding 
an image at successive intensity levels, and following the "parent"-"daughter" re- 
lationship between the structures identified at the different levels. Clearly, the two 
techiques applied to the same data produce very different sets of clumps. 

This variety of procedures for defining clumps illustrates the ambiguity inherent 
in defining a finite object that is actually part of a continuum, and implies that the 
very concept of a clump's mass carries with it a certain level of inherent uncertainty. 



Clump masses evolve in time 

A turbulent density fluctuation is a local density enhancement produced by a ve- 
locity field that at some moment in time is locally convergent, as indicated by eq. 
( [Tol l. This process accumulates mass in a certain region of space ("the clump"), with 
the natural consequence that the clump's mass must increase with time, at least ini- 
tially, if the clump is defined, for example, as a connected object with density above 
a certain threshold. This definition corresponds, for example, to clumps defined as 
compact objects observed in a particular molecular tracer, since such tracers require 
the density to be above a certain threshold to be excited. 

The growth of the clump's density (and mass) lasts as long as the total pressure 
within the clump (which may include thermal, turbulent and magnetic components) 
is smaller than the compression's ram pressure. However, in a single-phase medium, 
once the turbulent compression subsides, the clump, which is at higher density than 
its surroundings, and therefore also at a higher pressure, must therefore begin to 
re-exp and, unless it manages to become gravitationally unsta ble and proceed to col- 



lapse (IVazquez-Semadeni et al.L l2005at iGomez et aU 120071 see also ^3-3b . There 



fore, the mass above the clump-defining density threshold may begin to decrease 
again. (Again, see ! j4.2.1l for the case of clumps forming in multi-phase media.) As 
we shall see in ^3.4.31 the fact that clumps' masses evolve in time has direct im- 
plications for the amount of magnetic support that the clump may have against its 
self-gravity. 



3.4 Production of Density Fluctuations. The Magnetic Case. 

In the magnetized case, the problem of density fluctuation production becomes more 
complex, as the turbulent velocity field also produces magnetic field fluctuations. In 
this section, we discuss two problems of interest in relation to SF: The correlation 
of the density and magnetic fluctuations, and the effect of the magnetic field on the 
PDF of density fluctuations. 
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3.4.1 Density-Magnetic Field Correlation 



This is a highly relevant issue in relation to SF, as the "standard" model of magnetica lly 



regul ated SF (SMMRSF; see, e.g., the reviews bv lShu et al.L Il987t iMouschovias , 



1991b predicted that magnetic fields should provide support for the density fluctu- 
ations ("clumps") against their self-gravity, preventing collapse, except for the ma- 
terial that, through AD, managed to lose its support (see the discussion in fj5A2}. 
Thus, the strength of the magnetic field induced in the turbulent density fluctuations 
is an important quantity to determine. 

Under perfect field-freezing conditions, the simplest scenario of a fixed-mass 
clump threaded by an initially uniform magnetic field, and undergoing an i sotropic 
gravit ational contraction implies that the field should scale as B « p 2 / 3 ( MestelL 
1966b . since the density scales as p « R 3 , where R is the clump's radius, while 
the flux-freezing condition implies that B °< R~ 2 . The assumption that the clump is 
instead oblate, or disk-like, with the magnetic field providing support in the radial 
direction and thermal pressure provid ing support in the directio n perpendicular to 
its plane, gives the scaling B °^ p~ 1 / 2 ( Mouschovias , 19761 1991 ) 

In a turbulent flow, however, the situation becomes more complicated. "Clumps" 
are not fixed-mass entities, but rather part of a continuum that possesses random, 
chaotic motions. In principle, unless the magnetic energy is much larger than the tur- 
bulent kinetic energy, the compressive motions that form a clump can have any ori- 
entation with respect to the local magnetic field lines, and thus the resulting density 
enhancement may or may not be ac companied by a corresponding mag netic field 
enhancement (cf. ^3. lb . In particular. lPassot & Vazquez-Semadenil (120031 hereafter. 
PV03) studied this problem analytically in the isother mal case, by decomposin g 
the flow into nonlinear, so-called "simple" waves (e.g., Landau & Lifshitz , 1959b. 
which are the nonlinear extensions of the well known linear MHD waves (e.g^ Shu , 
1992), having the same three well-known modes: fast, slow, and Alfven (IMann , 
19951) . For illustrative purposes, note that compressions along the magnetic field 
lines are one instance of the slow mode, while compressions perpendicular to the 
field lines (i.e., magnetosonic waves) are an instance of the fast mode. As is well 
known, Alfven waves are transverse waves propagating along field lines, and carry 
no density enh ancement in th e linear regime, but this restriction does not hold in the 
nonlinear case (lMannlll995l) . 

For simplicity and insight, PV03 considered the so-called "1+2/3 -dimensional" 
case, also known as "slab geometry", meaning that all three components of vector 
quantities are considered, but considering their variation with respect to only one 
spatial dimension. They concluded that each of the modes is characterized by a 
different scaling between the magnetic pressure (« B 2 ) and the density, as follows: 



B 2 oc Cl - pp slow, 



B 1 
B 1 



fast, 
Alfven, 



(12) 
(13) 
(14) 
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where c\ is a constant, and y m is a parameter that ca n take values in the range (1/2,2) 
depending on the Alfvenic Mach number (see also lMcKee & Zweibel , 1995b. N ote 
that eq. ( [T2l i implies that for p > ci//3 the slow mode disappears ( Mannl 1995 ). so 
that only the fast and Alfven modes remain. Conversely, note that, at low density, the 
magnetic pressure due to the fast and Alfven modes becomes negligible in compar- 
ison with that due to the slow mode, which approaches a constant. This implies that 
a log-log plot of B vs. p will exhibit an essentially constant value of B at very small 
values of the density. In other words, at low values of the density, the domination 
of the slow mode implies that the magnetic field exhibits essentially no correlation 
with the density. 

PV03 were able to test these results numerically by taking advantage of the slab 
geometry, which allowed to set up waves propagating at well-defined angles with 
respect to the mean magnetic field, and therefore being able to isolate, or nearly iso- 
late, the three different wave modes. The left panel of Fig.|5]shows the distribution 
of points in the B 2 -s space for a simulation dominated by the slow mode, exhibiting 
the behavior outlined above. In contrast, the middle panel of Fig. [5] shows the dis- 
tribution of points in the same space for a simulation dominated by the fast mode. 




Fig. 5 Two-dimensional histograms of the grid cells in numerical simulations in the InB 2 (vertical 
axis)-lnp (horizontal axis) space. The values of the variables are in non-dimensional code units. 
Left: A slab-geometry numerical simulation by PV03 dominated by the slow mode, exhibiting the 
behavior indicated by eq. M2\ . Middle: A numerical simulation in slab geometry, also by PV03, 
dominated by the fast mode, exhibiting the behavior indicated by eq. dl3t . Right: A fully turbulent, 
therm ally-bistable, three-dimensional simulation of molecular cloud formation by Baneriee et al. 
l l2009h , in which all wave modes are active. 



The most important conclusion from eqs. <TT~2b — (TT~4T> is that each mode produces a 
different dependence of the magnetic field strength with the density. Therefore, in a 
turbulent flow in which all modes are active, the instantaneous value of the density 
at a certain location in physical space is not enough to determine the value of the 
magnetic field strength there. Instead, this value depends on the history of modes of 
the nonlinear waves that have passed through that location, naturally implying that, 
within a large cloud, a large scatter in the measured values of the magnetic field is 
expected, as illustrated in the right panel of Fig. [5] 
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Finally, at high densities, where the slow mode disappears and the fast and Alfven 
modes dominate and in general produce a magnetic pressure that does increase with 
increasing density, one expects a better-defined correlation between density and the 
magnetic field, although with significant scatter, since the dependence of B 2 with p is 
different for the two modes. These results are in qualitative agreeme nt with detailed 



statistical analyses of the magnetic field distribution in the ISM (Crutcher et al 
201(1 



3.4.2 Effect of the Magnetic Field on the Density PDF 



According to the discussion in £|3 .3-11 the dependence of the pressure on density 
determines the shape of the density PDF, being a lognormal for the isothermal case, 
7c = 1 ■ In the presence of a magnetic field, it would be natural to expect that the mag- 
netic pressure, which in general does not need behave as an isothermal polytrope, 
might cause deviations from the lognormal density PDF associated to isothermal 
turbulent flows. 

However, the discussion above on the density-magnetic field correlation (or 
rather, lack thereof), implies that the magnetic pressure does not have a system- 
atic effect on density fluctuations of a given amplitude, as the value of the mag- 
netic field is not uniquely determined by the local value of the density. PV03 con- 
cluded that the effect of the magnetic pressure was more akin to a random forc- 
ing in the turbulent flow than to a systematic pressure that opposes compression. 
As a consequence, the underlying density PDF determined by the functional form 
of the thermal pressure did not appear to be significantly affected by the presence 
of the magnetic field, except under very special geometrical setups in slab geom- 
etry, and that in fact are unlikely to persist in a more general three-dimensional 
setup. The persistence of the underlying PDF dictated by the thermodynamics in 
the presence of the magnetic field is in agreement with numerical studies of isother- 
mal MHD turbulent flows that indeed have found approximately lognorma l PDFs 



mal MHD turbulent flows that indeed nave round app : 
in isothermal flows (e.g., Padoan & Nordlundi 19991 : 



Ostriker et all 119991 12001 . 

Vazquez-Semadeni & GarcfaT 200 it Beresnvak et al. , 2005b . and bimodal density 



PDFs in thermally-bistable flows dGazol et all 20091) . which we discuss in ^4. 1 



3.4.3 Evolution of the Mass-to- Magnetic Flux Ratio 



The discussion in i j3.3.2| implies that the mass deposited in a clump by turbulent 
compressions is a somewhat ill-defined quantity, depending on where and how one 
chooses to define the clump's "boundaries", and moreover is time-dependent. This 
has important implications for the so-called mass-to-flux ratio (MFR) of the clump, 
and therefore, for the ability of the magnetic field to support the clump against its 
self-gravity. 



As is well known (iMestel & Spitzen.119561) . a virial balance analysis implies that, 
for a cloud of mass M threaded by a uniform field B, gravitational collapse can only 



16 



Enrique Vazquez-Semadeni 



occur if its MFR satisfies 



cm 



where a is a constant of order unity w hose precise value depends on the shape and 



mass distribution in the cloud (see, e.g jMestel & Spitzeiill956tlNakano & Nakamura , 



1978c IShul 119921) . Otherwise, the cloud is absolutely supported by the magnetic 
field, meaning that the support holds irrespective of the cloud's density. In what fol- 
lows, we shall denote the MFR, normalized to this critical value, by jU. Regions with 
jU > 1 are called magnetically supercritical, while those with ji < 1 are termed mag- 
netically subcritical. Traditionally, it has been assumed that the cloud's mass is well 
defined. However, our discussions above ( ^3.3.25 suggest that it may be convenient 
to revisit these notions. 

When considering density enhancements ("cores") formed by turbulent compres- 
sions within a cloud of size L, it is convenient to assume that the initial condition 
for the cloud is one with uniform density and magnetic field, and that the turbulence 
produces local fluctuat ions in the density and the magne tic field strength. A simple 



argument advanced by IVazquez-Semadeni et al.l (12005 al) then shows that the MFR 



of a core of size I and MFR )X{, must be within the range 

£ 

IM>-<He< to, (16) 

where fiQ is the MFR of the whole cloud. The lower limit applies for the case when 
the "core" is actually simply a subregion of the whole cloud of size L, with the 
same density and magnetic field strength. Since the density and field strength are 
the same, the core's mass simply scales as (£/L) 3 , while the magnetic flux scales 
as (£/L) 2 . Therefore, the MFR of a subregion of size £ scales as (£/L). Of course, 
this lower-limit extreme, corresponding to the case of a "core" of the same density 
and field strength as the whole cloud, is an idealization, since observationally such 
a structure cannot be distinguished from its parent cloud. Nevertheless, as soon as 
some compression has taken place, the core will be observationally distinguishable 
from the cloud (for example, by using a tracer that us only excited at the core's 
density), and the measurement of the MFR in the core will be bounded from below 
by this limit. 

On the other hand, the upper limit corresponds simply to the case where the entire 
cloud of size L has been compressed isotropically to a size £, since in this case both 
the mass and the magnetic flux are conserved, and so is the MFR. 

This reasoning has the implication that the MFR that is measured in a core within 
a cloud must be smaller than that measured for the whole cloud, as long as the con- 
dition of flux-freezing holds. Note that one could argue that this is only an observa- 
tional artifact, and that the physically relevant mass is that associated to the whole 
flux tube the core belongs to, but this is only reflecting the ambiguity discussed 
above concerning the masses of clumps. In practice, the physically relevant mass 
for the computation of the MFR is the one responsible for the local gravitational 
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potential well against which the magnetic field is providing the support, and that is 
precisely the core's mass only, not the mass along the entire flux tube, especially 
when phase transitions are involved. 
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Fig. 6 Evolution of the mass-to-flux ratio (MFR, /I), normalized to the critical value, in 
cores formed in numerical simu lations of continuously driven MHD isothermal turbulence by 
IVazque z-Semadeni et al. ( 20053). The cores are defined as connected sets of grid cells with density 
above a density threshold « t , and are followed over time. Three values of n t are used, illustrating 
how the value of the MFR depends and evolves as successively more internal regions of the density 
fluctuation are considered. Left: Evolution of /i in a core that does not collapse. At all values of n t , 
the MFR first increases and then decreases again. Right: A core that does collapse. In this panel, 
two lines are shown at the largest value of « t , because the parent clump splits into two cores at 
this threshold. In both panels, the rate of variation of /I is larger for larger values of nt, and jl for 
the innermost region («, = lOO/io, where «rj is the simulation's mean density) is seen to start lower 
than that of the envelope, and to overtake it as the degree of mass concentration is increased, in this 
case as a consequence of numerical diffusion, which plays a role analogue to that of AD. 



Note also that, as discussed in £|3 . 3 .21 the clump's mass must be evolving in time. 
If the compression is occurring mostly along field lines (since in this direction the 
magnetic field presents no resistance to it), then the magnetic flux remains roughly 
constant, while the mass increases at first, and later it possibly decreases if the clump 
begins to re-expand. Otherwise, if the core becomes massive enough that it becomes 
gravitationally unstable and supercritical, it must begin to collapse gravitationally 
(Fig. [6]). At this point, the rapid density enhancement at the core in turn enhances 
the action of AD (e.g. JShu et al.Ul987tlMouschoviasl.ll99l[ see also £13.11 ). causing 
the magnetic flux to escape the core0 so that the latter eventually acquires a larger 
value of /i than its envelope. Thus, the prediction from this dynamic scenario of core 



3 Note that this "escape" is meant in a Lagrangian sense, i.e., following the flow. That is, con- 
sidering a certain fluid parcel as it contracts, AD causes the flux to be "left behind" from the fluid 
particles that make up the parcel. Conversely, in an Eulerian sense, the magnetic flux remains fixed, 
but the fluid parcel increases its mass in this frame. 
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formation is that cores in early stages of evolution should exhibit smaller values of 
the MFR than their envelopes, while cores at more advanced stages should exhibit 
larger values of the MFR than t heir envelopes. Eviden ce in this direction has begun 



to be collected observationally (ICrutcher et all I200S ) , as w ell as through synthetic 



observations of numerical simulations dLunttila et al 



2009). 



Before closing this section, an important remark is in order. Recent numerical and 
observational evidence (cf. Sec I5.41 i suggests that the "turbulence" in star-forming 
molecular clouds may actually consist of a hierarchy of gravitational contraction 
motions, rather than of random, isotropic turbulence. In such a case, the physical 
processes discussed in this section are still applicable, with the precision that the 
converging flows that produce the clumps may be driven by larger-scale gravita- 
tional collapse rather than by random turbulent compressions, and the only part of 
the previous discussion that ceases to be applicable is the possibility that some cores 
may fail to collapse and instead re-expand. If the motions all have a gravitational ori- 
gin, then all cores must be in their rout to collapse. It is worth pointing out that in 
this case, what drives the collapse of an apparently subcritical core is the collapse 
of its parent, supercritical structure. 



4 Turbulence in the Multiphase ISM 



In the previous sections we have separately discussed two different kinds of physi- 
cal processes operating in the ISM: radiative heating and cooling, and compressible 
MHD turbulence in the special case of isothermality. However, since both operate 
simultaneously, it is important to understand how they interact with each other, es- 
pecially because the mean density of the Galactic ISM at the Solar galactocentric ra- 
dius, (n) ~ 1 cm~ 3 , falls precisely in the thermally uns table range. This problem has 
been investigated numerically by various groups (e.g., Hennebelle & Perauhj 19991 



been investigated numerically by various groups (e.g., Hennebelle & Ferault, lyyy, 
2000tlWalder & Folini|.l200oTlKovama & Inutsukall200dl2002tlVazquez-Semadeni et al 



2000ai l2003ai 1200611200711201 it iHennebelle etall 120081: iBaneriee et all 12009 1 : 



Gazol et all |200l[ l2005t |2009t iKritsuk & Normani 12002 



2002 : Piontek & Ostrikerl 



2004, 2005; Audit & Hennebelle, 



Sanchez-Salcedo et al 
20051120 lOtlHeitsch et al 



2005ilHennebelle & AuditL l2007l)T and in this section we review their main results. 



4.1 Density PDF in the Multiphase ISM 

A key parameter controlling the interaction between turbulence and net cooling is 
the ratio r\ = T c / T t , where T c rj e / (ji m mnpA ) is the cooling time and T t w L/U is 
the turbulent crossing time. The remaining symbols have been defined above. In the 
limit tj ^ 1, the turbulent compressions' dynamical evolution occurs much more 
rapidly than they can cool, and therefore the compressions behave nearly adiabati- 
cally. Conversely, in the limit 77 <C 1, the fluctuations cool down essentially instan- 
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taneously while the turbulent compression is evolving, and thus they t end to reach 
the thermal equilibrium pressure P eq as soon as they are producecfl ( Elmegreeri 



19911 Passot et al.L Il995t ISanchez-Salcedo et al.L 120021 ; IVazquez-Semadeni et al 



2003at iGazol et aU 12005b . Because in a turbulent flow velocity fluctuations of a 
wide range of amplitudes and size scales are present, the resulting density fluctua- 
tions in general span the whole range between those limits, and the actual thermal 
pressure of a fluid parcel is not uniquely determined by its density, but rather de- 
pends on the details of the velocity fluctuation that produced it. This causes a scatter 
in the values of the pressure around the thermal-equilibrium value in the pressure- 
density diagram (Fig. [7] left panel), and also produces significant amounts of gas (up 



forbidden thermally unstable range (Gazol et all 2001 ; de Avillez & Breitschwerdt, 


2005; Audit & Hennebelle, 2005 


; MacLowetal.. 2005 


), a result that has been en- 


countered by various observational studies as well (e.g., 


Dickevetal.. 1978; Heiles, 



20011) . In any case, the tendency of the gas to settle in the stable phases still shows 
up as a multimodality of the density PDF, which becomes less pronounced as the 
rms turbulent velocity increases (Fig. [7] right panel). 




Fig. 7 Left: Two-dimensional histogram of the grid cells in the pressure-density diagram for a two- 
dimensional simulation of turbulence in the thermally-bistable atomic medium, with rms velocity 
dispersion of 9 km s , a numerical box size of 100 pc, and the turbulent driving applied at a scale 
of 50 pc. Right: Density PDF in simulations like the one on the left panel, but with three different 
values of the rms velocity: 4.5 km s (solid line), 9 km s (dotted line), and 1 1.3 km s (dashed 
line). The peaks in the P DF are seen to to become less pronounced as the rms velocity increases. 
From lGazol et alj j2005l) . 



It is often believed that fast cooling directly implies isofhermality. However, this is a misconcep- 
tion. While it is true that fast cooling is a necessary condition for approximately isothermal be- 
havior, the reverse implication does not hold. Fast cooling only implies an approach to the thermal 
equilibrium condition, but this need not be isothermal. The precise form of the effective equation 
of state depends on the details of the functional dependence of the heating and cooling functions 
on the density and temperature. 
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4.2 The Formation of Dense, Cold Clouds and Clumps 
4.2.1 The Non-Magnetic Case 



A very important consequence of the interaction of turbulence (or, more generally, 
large-scale coherent motions o f any kind) and TI is that the former may nonlin- 
early induce the latter. Indeed, iHennebelle & Peraulj (119991 see also Koyama & 
Inutsuka 2000) showed that transonic (i.e., with M s ~ 1) compressions in the WNM 
can compress the medium and bring it sufficiently far from thermal equilibrium that 
it can then undergo a phase transition to the CNM (Fig. [8] left panel). This process 
amounts then to producing a cloud with a density up to lOOx larger than that of the 
WNM by means of only moderate, transonic compressions. This is in stark contrast 
with the process of producing density fluctuations by pure supersonic compressions 
in, say, an isothermal medium, in which such density contrasts would require Mach 
numbers M s ~ 10. It is worth noting t hat the turbulent velocity dispersion of ~ 8- 



1 1 km s 1 in the warm Galactic ISM dKulkarni & HeilesL Il987t iHeiles & Troland , 



2003) is, precisely, transonic. 
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Fig. 8 Left: Evolutionary path (dashed line) in the P vs. p diagram of a fluid parcel initially in the 
WNM after suffering a transonic compression that nonlinearly triggers TI. The solid and dotted 
lines show the locus of P eq (p), the solid sections corresponding to linear stability and the dotted 
ones to linear instability. The solid section to the left of the dotted line corresponds to the WNM 
and the one at the right, to the CNM. The pe rturbed parcel evolves from left to right along the 
dashed line. From [Koyama & Inutsukal {2000,). Right: Projected (or column) density structure of 
the resulting GMC in a numerical simulation of its formation by colliding WNM streams, and its 
subsequent evolution. The numerical box has a size of 15 pc on a side, and the resolution is 1200 3 
grid cells. The "GMC" is seen to consist of the agglomeration of a huge number of small clumps, 
which have fragmen t ed du e to the action of combined instabilities in the compressed gas. From 
lAudit & Hennebelld l2O10h . 



Moreover, the cold clouds formed by this mechanism have typical sizes given by 
the size scale of the compressive wave in the transverse direction to the compres- 
sion, rather than having to be of the same size scale as the fastest growing mode of 
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TI, which is very small (~ 0. 1 pc; cf. jQ. The ini t ial sta ges of this process may pro- 
duce thin CNM sheets (jVazquez-Semadeni et al. . 120061) . which are in fact observed 
dHeiles & Trolandl, 120031) . However, such sheet s are quickly dest abilized, apparently 



by a combination of nonlinear thin shell (NTSI: I Vishniac L 1 1 9941) . Kelvin-Helmholtz 



and Rayleigh-Taylor instabilities dHeitsch et all 120051) . fragmenting and becoming 



al.. 



turbulent. This causes the clouds to become a complex mixture of cold and warm 
gas, where the cold gas is distributed in an intrincate network of sheets, filaments 
and clumps, possibly permeated by a dilute, warm background. An example of this 
kind of structure is shown in the right panel of Fig. [8] 

A noteworthy feature of the clouds and clumps formed by TI is that, contrary to 
the case of density fluctuations in single-phase media, they can have more clearly 
defined and long-lasting boundaries. This is because their boundaries may be de- 
fined by the locus of the interface between the cold and warm phases, which, once 
formed, tends to persist over long timescales compared to the dynamical time, be- 
cause the two phases are essentially at the same pressure. Under quasi-hydrostatic 
conditions, these boundaries hav e little or no mass flux across them (i.e., they are 



contact discontinuities; see, e.g., IShuL 119921) . with whatever little mass exchange 



that happens being due to evaporation or condensation, when the thermal pressure 
differs from the saturation value between the phases (e.g..lZerD ovich & Pike FNei , 
19691; IPenston & Brownlll970tlNagashima et all [2005t llnoue et al.ll2006l) . The lat- 



ter two papers have in fact proposed that such evaporation may contribute to the 
driving of interstellar turbulence, although the characteristic velocities they obtained 
( 1 km s~ 1 ) appear to be too small for this to be the dominant mechanism for driv- 
ing the large scale ISM turbulence, with characteristic speeds of ~ 10 km s _1 . 

However, in the presence of large-scale (> 10 pc) and large-amplitude ( £ 10 km 
s ) motions, corresponding either to the supernova-driven global ISM turbulence , 



to larger-scale instabili ties, such as the magneto- Jeans (e.g.. lKim & Ostrikeii 120011) 
or magneto-rotational ( Balbus & Hawlevl 1991 ) ones, or simply to the passage of 
sparal arms, the nonlinear triggering of TI implies that the phase transition fronts 
bounding the clouds and clumps are not contact discontinuities, but rather phase 
transition fr onts - structures analogous to shocks, through which a substantial mass 
flux occurs dVazquez-Semadeni et al. . 120061: iBanerjee et al.l . 120091) . This is mecha- 
nism is illustrated in Fig. [9] which clearly shows the rapid growth of a clump in a 
numerical simulation by accretion of diffuse material. This is in stark contrast with 
classical i deas that the clumps grew by coagulation on very long timescales (~ 100 
Myr; e.g., lKwanl[i979l) . 

It should be emphasized, however, that the density does not necessarily always 
present a jump at the clump boundaries. As discussed in iH.ll the presence of tur- 
bulence in the diffuse medium also implies a certain degree of mixing, and the ex- 
istence of a certain fraction of the mass that is traversing the unstable range. In Fig. 
|9]this can be observed as the greenish regions, especially in the left panel. 

We conclude then that, although in thermally bistable flows clump boundaries are 
in general better defined than in turbulent isothermal flows because of the density 
jumps induced by the thermal bistability, this does not imply that they are impene- 
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Fig. 9 Density cuts through the plane x = 2. 5 pc in an adaptive-mes h refinement (AMR) numerical 
simulation of molecular cloud formation by iBaneriee et alj j2009h . illustrating the clump growth 
mechanism. The numerical box size is 256 pc, and the maximum resolution is 0.06 pc. The arrows 
show the projection of the velocity field on this plane. Left: The clump at time t = 21.6 Myr. Right: 
The clump at t = 22.5 Myr. Comparing the two times the growth of the clump is evident. The 
velocity field is seen to generally point towards the clump, indicating that material from the diffuse 
external medium is entering the clump, causeing its growth. 



trable boundaries that restrict the flow of the medium. Rather, the clumps are formed 
and then grow by accretion of diffuse material across these phase transition fronts. 

Finally, it is important to n ote that , contrary to earlier ideas of GMC formation 
by coagulation of clumps (e.g. iKwanl 1 1 979t blitz & Shul [l98oh . in the scenario of 
GMC formation described above, the compressions in the WNM tend to initially 
form thin sheets of CNM, w hich then fragment and thicken due to the combine d 
action of various instabilities dVazquez-Semadeni et al. ■ l2006tlHeitsch et alll2006h . 
so that GMCs may actually consist of huge conglomerates of small clumps, as il- 
lustrated in the rig ht panel of F ig. [8] This is consistent with the observed clumpy 



structure of GMCs (Blitz, 1993, sec. VII) 



4.2.2 The Magnetic Case 

In the presence of a magnetic field, the process of cloud formation by phase 
transitions to the cold phase requires further considerations. First, the orientation 
of the compressive motion relative to that of the magnetic field strongly influ- 
ences the ability of the com pression to trigger a transition to the dense phase. 



Hennebell e & Perault d2000h investigated this problem by means of two-dimensional 



simulations, finding that, for a certain value of the magnetic field strength, and a 
given sonic Mach number of the compression, there exists a maximal angle be- 
tween the direction of compression and the direction of the magnetic field beyond 
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which no phase transition is induced. They found this angle to typically lie between 
20 to 40 degrees, for typical val ues of the warm neutral medium. 



Hennebell e & Perault also found that, when the formation of a cloud does 



occur, either the field is re-oriented along the compression (in the case of weak 
fields), or the flow is re-oriented along field lines (in the case of strongter fields), and 
the accumulation of gas to form the clump ends up being aligned with the magnetic 
field. In addition, Inoue & Inutsukal d2008l) have found that compressions perpen- 
dicular to the magnetic field strongly inhibit the formation of dense, molecular-type 
clouds, and that, in this case, only diffuse HI clouds manage to form. As a con- 
sequence, the discussion of cloud formation can be made in terms of compressions 
parallel to the magnetic field without loss of generality. We will take up this problem 
again in i)5.4| when we discuss the onset of gravitational collapse of the clouds. 



5 The Nature of the Turbulence in the various ISM Components 



5.1 Generalities 



As discussed in the previous sections, the ionized and atomic components of the 
ISM consist of gas in a wide range of temperatures, from T ~ 1 6 K for the HIM, to 
T ~ 40 K for the CNM. In particular. iHeiles & Trolandl (120031) report temperatures 



in the range 500 < T < 10 4 K for the WNM, and in the range 10 < T < 200 K for 
the CNM. The WIM is exp ected to have T ~ 10 4 K. T his implies that the adiabatic 
sound speed, given by (e.g jLandau & Lifshitd. [l959b 



/ ykT 



10.4 km s" 



(io 4 k) 



1/2 



(17) 



will also exhibit large fluctuations in the medium. In the following sections we dis- 
cuss the implications of these ranges for the various ISM components. 



5.2 The Warm Ionized Medium 



Collecting measurements of interstellar scintillation (fluctuations in amplitude and 
phase of radi o waves caused by sca ttering in the ionized ISM) from a variety of 
observations. lArmstrong et al.l ( 11995b estimated the pow er spectrum of densit y fluc- 
tuations in the WIM, finding that it is consistent with a Kolmogorov jl94ll) spec- 



trum, a result expected for weakly compressible flows dBavlv et al 
10 8 < L < 10 15 



1992), on scales 



cm. 



Mo re recently, using data from the Wisconsin Ha Mapper Observatory, Chepurnov & Lazarianl 
(2010) lave been able to extend the spectrum to scales ~ 10 19 cm, suggesting that 
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the WIM behaves as an incompressible turbulent flow over size scales spanning 
m ore than 10 order s of magnitude. This suggestion is supported also by the results 
of lHill et al.1 (12008b who, by measuring the distribution of Ha emission measures in 
the WIM, and comparing with numerical simulations of turbulence at various Mach 
numbers, concluded that the sonic Mach number of the WIM should be ~ 1.4-2.4. 
Although the WIM is ionized, and thus should be strongly coupled to the mag- 
netic field, the turbulence then being magnetohydrodynamic (MHD), Kolmogorov 
scaling should still apply, accordin g to the theory of incompressible MHD fluctu- 
ations dGoldreich & Sridhari 1 19951) . The likely sources of kinetic energy for these 
turbulent motions are stellar energy sources such as supernova explosions (see, e.g., 



Mac Low & Klessen, 2004). 



5.3 The Atomic Medium 

In contrast to the relatively clear-cut situation for the ionized ISM, the turbulence 
in the neutral (atomic and molecular) gas is more complicated. According to the 
discussion in $5.1\ the temperatures in the atomic gas ma y span a continuous range 



from a few tens to several thousand degrees. Additionally. iHeiles & Trolandl d2003l) 
report column density-weighted rms velocity dispersions a v ~ 11 kms -1 for the 
WNM, and typical internal motions of M s ~ 3 for the CNM. It is thus clear that 
the warm, or diffuse, gas is transonic (M s ~ 1), while the cold, or dense, gas is 
moderately supersonic. This occurs because the atomic gas is thermally bistable, 
and because transonic compressions in the WNM can nonlinearly induce TI and 
thus a phase transition to the CNM ( 34. U . the neutral atomic medium is expected 
to consist of a complex mixture of gas spanning over two orders of magnitude in 
density and temperature. 



It is worth noting that early pressure-equilibrium models (e.g. jField et al.lll969 



McKee & Ostrikerl 1 1977ft proposed that the unstable phases were virtually nonexis- 
tent in the ISM, but the observational and numerical results reported in S|4T|suggest 
that a significant fraction of the atomic gas mass lies in the unstable range, as it is 
transiting between the stable phases. Also, numerical simulations of such systems 
suggest that the velocity dispersion within the densest clumps is subsonic, but that 
the velocity dispersion of the clumps within the diffuse substrate is supersonic with 
resp ect to the clumps' sound speed (although subson ic with respect to the warmest 



gas: lKovama & Inutsukal, 120021: iHeitsch et al.L 120051) . 
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5.4 The Molecular Gas 



5.4.1 Molecular Clouds: Supersonically Turbulent, or Collapsing? 



Mole cular clouds (MCs) have long been kn own to be stron gly self-gravitatin 
(e.gjGoldreich & Kwari 1 1974 lLarsonl[l98ll) . In view of this, iGoldreich & Kv ... 
(1974) initially proposed that MCs should be in a state of gravitational collapse, 
and that the observed motions in MCs (as derived by the non-thermal linewidths 
of molecular lines) corresp onded to this collapse. However, shortly thereafter, 



Zuckerman & Palmer! (119741) argued against this possibility by noting that, if all the 
molecular gas in the Galaxy, with mean density n ~ 100 cm~ 3 and total mass M mo \ ~ 
were in free-fall, then a simple estimate of the Galaxy's SF rate (SFR), given 
Mnoi/% ~ 200 Mr.) yr~\ where % = \/ 3n/32Gp is the free-fall ti me, 



by SFR 



would exceed the observed rate of ~ 2 M© yr _1 (e.g., Chomiuk & PovicR 201 lb by 
about two orders of magnitude. This prompted the suggestion (iZuckerman & Evans , 
19741) that the non-thermal motions in MCs corresponded instead to small-scale (in 
comparison to the clouds' sizes) random turbulent motions. 

The need for these motions to be confined to small scales arose from the need 
of the turbulent (ram) pressure to provide an approximately isotropic pressure that 
could counteract the clouds' self-gravity at large, thus keeping them from collaps- 



ing and maintaining them in near virial equilibrium ( Larsonl 



1981). Because tur- 



Landau & Lifshitzll 1959b 



bulence is known to be a dissipative phenomenon (e.g 
research then focused on finding suitable sources for driving the turbulence and 
avoiding rapid dissipation. The main driving source was consid e red to be energy 



inject i on from stars (e . g., Norman & Silki 1980b McKeel Il989t iLi & Nakamural 

al. I 120061 : ~ 



2006t iNakamura & LiL 120071 iKrumholz et al 



2006; Carroll et al, 2009, 2010; 



Wang et alT 20 id see also the reviews by Mac Low & Klessen 2004 and Vazquez 
Semadeni 2010), and reduction of dissipation was proposed to be accomplished by 
having the turbulence being MHD, and consist ing mostly of Al fven waves, which 



were thought not to dissipate as rapidly (e.g.. IShu et al 
provide an isotropic pressure ( McKee & ZweibeT 1995 ). 



1987), and which could 



However, in the last decade several results have challenged the turbulent pressure- 
support scenario: 1) Turbulence is known to be characterized by having the largest- 
velocities occur at the largest scales, and MCs are no exception, exhibiting scal- 
ing relations between velocity dispersion an d size suggest i ng that the largest ve- 
locity differences occur at the largest scales ( Larsonl 1981 : Heyer & BruntL 2004 



Brun t et al. . I2009L Fig. \W[ left and middle panels). This is inconsistent with the 



small-scale requirement for turbulent support. 2) It was shown by sever al groups that 



MHD turbulence dissipates just as rapidly as hydrody namic turbulence ( IMac Low et al 



1998; IStone et all 119981; iPadoan & Nordlu nd. 1999), dismissing the notion of re 



duced dissipation in "Alfven-wave turbulence", and thus making the presence of 
strong driving sources for the turbulence an absolute necessity. 3) Clouds with 
very different contributions from various turbulence-driving mechanisms, includ- 
ing those with little or no SF activity, such as the so-called Maddal ena's cloud 
(Maddal ena & Thaddeusll 1985b . show similar turbulence characteristics (Wi lliams et al 
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1994; lSchneider et aUl201 II) . suggesting that stellar energy injection may not be the 



main source of turbulence in MCs. 



t = 20.00 Myr 




Fig. 10 Left and middle panels: Second eigenimages obtained by Principal Component Analysis 
of spectroscopic data of the star-forming region NGC 1333, showi ng the main contribution to the 
linewidth of molecular emission in this region lBruntetalll2009l) . The middle image shows the 
region enclosed in the rectangle i n the left image. B lack and white colors represent oppositely- 
signed components of the velocity. Brunt et ah (2009) describe the pattern as a "dipole", in which 
large-scale patches of alternating velocity direction are observed. This is seen in both the large- 
scale (left) and the small-scale (middle) images. Right panel: Image of the projected density field 
of a 3D numerical simulation with cooling, self-gravity, and magnetic fields, representing the for- 
mation of a dense atomic cloud by the collision of WNM streams in the direction perpendicular to 
the plane of the figure. The time shown is 20 Myr after the start of the simulation. The black dots 
denote "sink" particles, which replace local collapsing zones in the simulation. The whole cloud is 
also collap sing, although its collapse is not completed yet by the end of the simulation, at / = 31 
Myr. From Vazquez-Semade ni et al.1 1201 ll) . 



Moreover, simulations of dense cloud formation in the nonmagnetic case have 
shown that, once a large cold CNM cloud forms out of a collision of WNM 
streams, it quickly acquires a large enoug h mass that it can begin to collapse 
gravitationally in spite of it being turbulen t (fV azauez-Semadeni et al. . 2007 . 20101 : 



Heitsch & Hartmanrl l2008t iHeitsch et all l2008a ). The enhance ment in its col 



umn d ensity promo t es the formation of molecular hydr ogen (H2) dHartmann et al 



200 U iBergin et all 120041 IHeitsch & Hartmannl 120081) . Thus, it appears that the 



formation of a molecula r cloud may require previous gravitational contraction 
(see also iMcKee , 119891) . In addition, according to the discussion in Sees. 14.2 
and 15.31 the CNM clouds formed by converging WNM flows should be born 
turbulent and clumpy. This turbule nt nature of the clo u ds fur ther promotes the 
formation of molecular hydrogen ( Glover & Mac Lowi 20071). T he simulations 
( Vazquez-Semadeni et al. , 2007 , 20 lot Heitsch & Hartmannl 2008 ) show that the 
nonlinear, turbulent density fluctuations can locally complete their collapse before 
the global collapse of the cloud is completed (Fig. [lOl right panel), both because 
their densities are large enough that their free-fall time is significantly shorter than 
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that of the whole cloud (IHeitsch et al. . l2008bt iPon et all 1201 lb . and because the 
free-fall time of a flattened or elongated cloud may be much larger th an that of an 



approximately isotrop ic clump of the same volume density within it (IToala et al 



2012tlPon et all 120 12b 



In addition, the turbulent velocities initially induced in the clouds by the con- 
verging flows in the simulations are observed to be relatively small (only mod- 
erately supersonic [M s ~ 3] with respect to the dense gas), while strongly su- 
personic (M s ^10) velocities like those observed in real molecular clouds only 
develop later, due to the ensuing gravitational contraction. Indeed, CNM clouds 
clouds are observed to typically h ave moderately supersonic (M s ~ 3) velocity 
dispersions dHeiles & Trolandil2003b . while GMCs are observed to have much larger 
turbulent rms Mach numbers, M s ~ 10- 20. Also, the clumps w ith highest internal 
velocity dispersions in the simulations of lBaneriee et al.l (12009b were those that had 
already formed collapsed objects ("sink" particles), even though energy feedback 
from the sinks was not included, implying again that the largest velocities develop 
by the action of self-gravity. 

Finally, it is very important to note that the possibility of MCs being in gravita- 
tional co llapse is not in contradiction with any observed properties of MCs since, as 
noted by Ballesteros-Paredes et al. ( 201 lab , the magnitudes of the virial and free-fall 
velocities for a self-gravitating object are observationally indistinguishable. Thus, 
the interpretation of cloud energetics in terms of virial equilibrium is completely in- 
terchangeable by an interpretation of collapse. The main arg ument against the pos- 
sibility of gravitational collapse of MCs is the argument by Zuckerman & Palmer 
(119741) that it would lead to exceedingly large SFRs. We discuss the possible reso- 
lution of this conundrum in £ 16.21 

All of the above evidence suggests that the observed supersonic motions in 
molecular clouds may have a significant, perhaps dominant, component of infalling 
motions, with a (probably subdominant) sup erposed random (turbulent) component 



remain ing from the initial stages of the cloud (IBate et all l2003: Ballesteros-Parede s et al 



201 lalbl), and perhaps somewhat amp lified by the collapse ( Vazquez-Semadeni et all 
19981; [Robertson & Goldreichi 1201 2) . In this scenario of hierarchical gravitational 



fragmentation, the first structures that complete their collapse are small-scale, high 
amplitude density fluctuations that are embedded withi n larger-scale ones, but of 
lower amplitude, which complete their collapse later (iVazquez-Semadeni et al 



2009b . The main role of the truly turbulent (i.e., fully random) motions is to provide 



the nonlinear density fluctuation seeds that will collapse locally once the global con- 
traction has caused their density to increase suff iciently for them to become locally 
gravitationally unstable ( Clark & BonnellL 12005 ). Evid ence for such multi-scale col 



lapse has recent l y begu n to be observationally detected (IGalvan-Madrid et all 12009 



Schneider etaD.l2010b 
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5.4.2 The Molecular Gas. Results Including the Magnetic Field 

According to the discussion in ^4.2.21 the formation of a cold dense atomic cloud 
involves accumulation of warm material along magnetic field lines in order for the 
necessary phase transition to the cold phase to occur. However, as discussed in the 
previous section, the formation of a molecular cloud probably requires the grav- 
itational contraction of the atomic cloud formed by the compression. Thus, in the 
presence of the magnetic field, this requires an understanding of the role of magnetic 
support; that is, of the evolution of the MFR. 

As is well known, and was reviewed in ^3.4.31 there exists a critical value of the 
MFR above which the magnetic field is no longer able to support the cloud against 
its own self-gravity. For accumulation along field lines, the criticality condition in 
terms of the mass column density E = pL and the field strength Bq for a cylindrical 



geometry is ( INakano & Nakamural.119781) . 



(Z/B ) crit « (4tt 2 G)- 1 / 2 « 0.16 G" 1 / 2 , (18) 

where p is the mass density and L is the cylinder length. This condition gives the 
accumulation length, i n terms of fiducial valu es representative of the ISM in the 



solar neighborhood, as (lHartmann et all 1200 ll) 



Bo \ 



-] 



l ^ 470 U^JIt^J pc ' (19) 

where we have assumed fi m = 1.27. In principle, if the Galactic field is primar- 
ily azimuthal, then the Galactic ISM at large is magnetically supercritical in gen- 
eral, because field lines circle around the entire Galactic disk, and thus sufficiently 
long distances are always available along them0Thus, the MFR of a system is not 
a uniquely defined, absolute parameter, but rather depends on where the system's 
boundaries are drawn. Also, recall that the critical value of the MFR depends on the 
local geometry of the system being considered. For instance, a system with spher- 
ical wmmetry has a critical value of (Z/Bo)crit = (67r 2 G)~ 1,/2 « 0.13 G -1 / 2 (e.g. 



Shull 19921) . somewhat smaller than that given by eq. ( flU . 

Now consider a cloud or clump that is formed by the accumulation of gas along 
field lines in general^ In the rest of this discussion, we will generically refer to the 
resulting density enhancement as a "cloud", referring to either a cloud, a clump, 
or a core. Although redistribution of matter along field lines does not in principle 
affect the total MFR along the full "length" of a flux tube, this length is a rather 
meaningless notion, since the flux tube may extend out to arbitrarily long distances. 



Note, however, that supercriticality does not necessarily imply collapse, since the gas may be 
thermally or otherwise supported, as is likely the case for the diffuse warm medium at scales of 
hundreds of parsecs. 

6 Since compressions perpendicular to the magnetic field cannot induce collapse of an initially 
subcritical region, as they do not change the MFR, and compressions o blique to the field can 
produ ce collapse by reorienting the directions of the flow and the field lines (Henn ebelle & Peraulj 
l200d) . our assumed configuration involves no loss of generality. 
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What is more meaningful is the MFR of the dense gas that makes up the cloud, 
since the cloud is denser than its surroundings, and thus it is the main source of the 
self-gravity that the field has to oppose. In fact, for the formation of a cloud out of 
flow collisi ons in the WNM, t he cloud's density is ~ 100 times larger than that of 



the WNM ( iField et al.L 119691: IWolfire et al.L 1 1995b . and so the latter's self-gravity 



is negligible. Thus, in this problem, natural boundaries for the cloud are provided 
by the locus of the phase transition front between the dense and the diffuse gas, 
allowing a clear working definition of the MFR. 

However, contrary to the very common assumption of a constant cloud mass, 
the formation of clouds by converging gas streams implies that the cloud's mass is 
a (generally increasing) function of t ime (cf. i ]3.3.21 i. a conclusion that has recently 
been reached observationally as well dFukui et al. 1 120091) . This means that, within the 
volume of the cloud, the MFR is also an increasing quantity, since the fl ux remains 



const ant if the flow is along field lines, while the mass increases (see also lShu et al 



20071) . If the cloud starts from essentially zero mass, this in turn implies that the 
MFR of a cloud is expected to start out strongly subcritical (when the cloud is only 
beginning to appear), and to evolve towards larger values at later times. Rewriting 
eq. (T% f or the column density, we see that the cloud becomes supercritical when 
(iVazquez-Semadeni et al.L 1201 11) 



1.5 x 10 



21 



Bo \ _ 
cm 



(20) 



where N = E/jimn is the number column density, and is to be measured along the 
field lines. The critical column density for magnetic criticality given by eq. (f20b 
turns out to be very similar, at least for solar neighbourhood conditions, to the crit- 
ical column density of hydrogen ato ms necessary for cold atomic gas to become 
molecular, Mr ~ l-2x 10 21 cm' 2 (e.g., Franco & Cox , 1986 ; van Dishoeck & BlackL 



1988tlvan Dishoeck & BlaQll998tlHartoann et aUl200ltlGlover & Mac Lowll2007llb 



Glover et al.l 1201(1 

Moreover, the critical column density given by eq. (|2Q] | is also very similar to 
that required for rendering cold gas gravitationally unstable, which is estimated to 
be 

n ( P/k ^ 1/2 



N, 



grav 



; 0.7 x 10 



(21) 



3000 K cm- 3 

dFranco & Cox , ^986t lHartmann et al U2001I) . Thus, the evolution of a cloud is such 

that it starts out as an atomic, unbound, and subrcritical diffuse cloud nVdz.quez.-Semadeni et al. 
200(\) and, as it continues to accrete mass from the warm atomic medium, it later be- 
comes mol ecular, supercritical, an d collapsing, all happening at approximately the 
same time ( Hartmann et al.ll200ll) . This is fully consistent with the observation that 
diffuse atomic clouds are in general strongly subcritical dHeiles & Trolandl 120051) 
and not strongly se lf-gravitating, while GMCs are approximately critical or moder - 
ately supercritical dCrutcherl.il 9991: iBourke et aUl200ltlTroland & Crutcheril2008l) . 
and are generally gravitationally bound (e.g., iBlitzL 1 1 993b . 
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It is important to note that this is in stark contrast to the SMMRSF, where it 
was considered that the magnetic criticality of a cloud was the main parameter de- 
termining whether it would form only low-mass stars and at a slow pace (in the 
case of subcritical clouds), or form clusters, including high-mass stars, and at a fast 
pace (supercritical clouds). This constituted a bimodal scenario of SF, and sub- and 
supercritical clouds constituted two separate classes. 

Instead, in the evolutionary scenario for MCs described above, clouds are ex- 
pected to evolve from being simultaneously atomic, subcritical and not strongly 
self-gravitating to being molecular, supercritical and strongly self-gravitating. Next, 
the roughly simultaneous transition to self-gravitating and supercritical, suggests 
that, in general, GMCs should be in a state of gravitational contraction, at least ini- 
tially, even in the presence of typical magnetic field strengths in the Galactic disk. 
Of course, significant scatter in the MFR is expected, bot h intrinsically (s ee ^3.4. \\ 
and as a consequence of observational uncertainties (e.g.. ICrutcheri 1 1 999b . and thus 
a certain fraction of the GMCs may remain subcritical up to significantly evolved 
stages, or even throughout their entire evolution. This case is discussed further be- 
low. 

The formation and evolution of mo lecular clouds in the magnetic ca se has been 



recently investigated numerically by iVazquez-Semadeni et al.l (1201 11) . who per 



formed simulations of GMC formation by compressions in the WNM aligned with 
the magnetic field, including ambipolar diffusion. Three cases were considered, one 
supercritical, with jU = 1.3, and two subcritical, with fx = 0.9 and 0,7, respectively. 
The initial magnetic field was considered uniform, and respectively having strengths 
of 2, 3, and 4 jj, G. In all cases, the mean density was 1 cm~ 3 and the temperature 
T = 5000 K. The compressions consisted of two oppositely-directed streams of gas 
at the mean density, and of length 1 12 pc, immersed in a 256-pc box. 

The evolution of the subcritical cases is worth discussing in detail, as it differs 
somewhat from simple expectations. These simulations produced a dense cloud that 
quickly began to contract gravitationally, similarly to non-magnetic simulations. 
This occurred because a uniform magnetic field does not provide any support, since 
the latter requires the existence of a magnetic gradient. Support builds up gradu- 
ally as the field lines are bent. The clouds thus contracted for a few tens of Myr, 
until the magnetic tension was large enough to halt the collapse, at which point 
they re-expanded, and entered an oscillatory regime, around the equilibrium con- 
figuration. However, due to the existence of diffusion (both numerical and from 
AD), local collapse eve nts managed to occur, in agreement with the notions from 
the SMMRSF (see also lMcKeel Il989l) . The notable difference, though, lied in the 
fact that the clouds only formed stars during the global contraction phase, especially 
at maximum compression, and essentially shut off in the re-expanding phase. This 
is in contrast to the SMMRSF, in which the GMCs at large were assumed to be in 
equilibrium and forming stars continuously, albeit slowly. 

The above discussion suggests that the possibility of star-forming molecular 
clouds being in a state of gravitational contraction, may hold even if they have sub- 
critical MFRs. The subsequent re-expansion of these clouds (or their remnants) may 
lead to a star-formation-inactive and quiescent phase, perhaps on their way to disper- 
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sal, if the clouds are exiting the spiral arms by that time, as proposed bv lElmegreen 
(120071) . 



6 Star Formation in the Turbulent ISM 

6.1 Does Molecular Cloud "Turbulence" Provide Support for 
Molecular Clouds? 



In the previous sections we have discussed how large-scale compressions in the gen- 
eral ISM produce density fluctuations, in particular by nonlinearly inducing phase 
transitions from the cold to the warm medium. Because the largest dimensions of the 
clouds thus formed are as large as the transverse dimension of the compression that 
formed them, they can soon find themselves being strongly gravitationally unsta- 
ble, and proceed to collapse. It is important to note that the large-scale compression 
forming the clouds may (and in fact, is likely to) have an origin different from the 
general turbulence in the ISM, such as, for exam ple, large-scale instabilities like the 
magneto-Jeans one (e.g. lKim & Ostrikeru200lh . or simply the passage of the stellar 
spiral-arm potential well. 

In fact, it is worth noting that just the turbulence driven by supernovae does 
not seem to be able to sustain itself, since the mass driven into a Jeans-unstable 
regime per unit time by the turb ulence is not enough to ma intain the same supernova 
rate that drives the turbulence ( Joung & Mac Low , 20061) . This conclusion is also 
supported by the fact that successive generations of triggere d SF do not appear to be 
able t o form stars as massive as in the previous generation dDeharveng & Zavagno , 
201 lb . Thus, it appears safe to conclude that the main driver of star formation is 
gravity at the largest scales. 

During the last decade, the main role of interstellar turbulence has been thought 
to be the regulation of SF, mainly on the basis of the assumption that the turbulent 
velocity dispersion contributes to the support of molecular clouds against their self- 
gravity, analogously to the role of th e thermal velocity dispersion, and perhaps in 



gravity, analogously to tne role ot tne thermal velocity dispersion, ana pernaps u 
eluding a scale-dependent amplitude dChandrasekharl 1951; Bonazzola etaUll98 



Bertoldi & McKeeLll992tlyazauez-Semadeni & GazolLll995l:lMac Low & Klesse 
2004|:lKrumholz & McKeell2005l:lHennebelle & Chabrierl.l2008l.l201 lUPadoan & Nordlundl 
201 lb . Thus, turbulence has been thought to provide support to clouds as a whole, 
while simultaneously inducing small-scale density fluctuations (clumps) within the 
clouds that may undergo gravitational collapse if they a re compressed enough for 



their Jeans mass to become smalle r than t heir actual mass (IVazquez-Semadeni et al 
2003bl iBallesteros-Paredes et all 120071) In this manner, global col lapse of the 



clouds could be prevented, avoiding the IZuckerman & Palmer! (11974) conundrum 
that the global collapse of molecular clouds would cause an excessively large SFR 
(see ^5.4. It . while at the same time allowing for the collapse of a small fraction 
of the mass, brought to instability by the local supersonic turbulent compressions. 
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However, this last notion was challenged by lOark & Bonnelll (120051) . who argued 
that the turbulence only provides the seed density fluctuations for subsequent grav- 
itational fragmentation, without significant local reductions in the Jeans mass in- 
duced by the turbulence. 



Moreover, iHeitsch & Hartmannl (|2008) showed that the fraction of mass with 
short free-fall times ( £ 1 Myr) in the clouds increases monotonically over time in 
the presence of self-gravity, indicating a secular evolution towards higher densities, 
while simulations with no self-gravity exhibited a stationary fraction of mass with 
short free-fall times, as would be the case in clouds supported against collapse by 
the turbulence (Fig.fTTTi. 
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Fig. 11 Time evolution of the mass (indicated by the color scale) at a given free-fall time (vertica l 
axis) in numerical simulations of cloud formation and evolution by Heitsch & Hartmann (2008). 
The top row shows the free-fall time in a linear scale, while the bottom row shows it in logarithmic 
scale. The panels on the left show a simulation with no self-gravity, while the middle and right 
panels show two different simulations with self-gravity. In the case with no self-gravity, the fraction 
of mass at a given free-fall time is seen to remain nearly constant, and the minimum free-fall time 
to remain at ~ 1 Myr, while in the cases with self-gravity, the minimum free-fall time decreases 
secularly. 



Finally, the simulations have also shown that the fraction of molecular gas 
also increases in time, so that the cloud would ind eed be classified as atomic 
in its early phase s, and as molecular in later ones ( Heitsch & Hartmannl 2008 



Clark et all 120121) . In particular, the latter authors have shown that the formation 
of CO-dominated regions only occurs ~ 2 Myr before SF starts, although signifi- 
cant amounts of H2 can appear earlier, although this component cannot be directly 
observed. 

All of the above evidence suggests that the strongly supersonic motions observed 
in MCs may be a manifestation of the gravitational contraction occurring in the 
clouds, rather than truly turbulent (i.e., random, roughly isotropic) motions, of a 
separate origin, that can counteract the clouds ' gravitational contraction. It is worth 
noting here that lKlessen & Hennebelle fcoich have recently shown that, in general, 
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the accretion power at scales from entire Galactic disks to protostellar disks, pass- 
ing through the GMC scale, is more than enough to drive the turbulence observed in 
these systems. However, it should be noted that, in the case of GMCs, this sugges- 
tion differs qualitatively from the nature of the motions discussed above. Rather than 
accretion driving turbulent motions in the clouds which can then support them, the 
discussion above suggests that the observed motions in the clouds are the infall it- 
self, with only a small, subdominant, truly random turbulent component superposed 
on them. In this case, these motions cannot provide support against the clouds' self- 
gravity. 



6.2 Regulation of Star Formation Via Feedback 



All of the above evidence strongly suggests that interstellar clouds undergo a secular 
evolution, starting their existence as moderatly supersonic, magnetically subcritical, 
sheet-like atomic clouds, and evolving towar ds becoming supercriti c al, mo lecular, 
collapsing objects. However, in this case, the IZuckerman & Palmer] (119741) SF co- 
nundrum (cf. jj5.4. It must be addressed. That is, if MCs are essentially in free-fall, 
how to prevent the SFR from being two orders of magnitude larger than it is ob- 
served to be in the Galaxy? 

Early studies proposed that ionizing radiation from massive stars should be able 
to disperse a cloud as early as whe n only ~ 10% of the cloud's mass has been con- 
verted to stars (see, e.g., sec. 4 of iFieldL [l970l) . so that the remaining 90% would 



be preven ted from forming an y more stars. This suggestion, however, was chal- 
lenged by Mouschovias ( 19761) . who argued that those estimates were based on the 
assumption of unrealistically low mean densities for the clouds (~ 10 cm -3 ), and 
that using more realistic values (~ 10 5 cm~ 3 ) would result in a grossly insufficient 
amount of ionization in the cloud, thus inva lidating the mechani sm as a suitable one 
for dispersing the clouds. As an alternative. iMouschoviasI (1 1 9761) proposed the basic 
notions for the SMMRSF: that the MCs should be magnetically subcritical in gen- 
eral, so that their envelopes would remain supported by the magnetic tension, while 
only the central core would be able to proceed to collapse through AD (cf. § 33.11 
and 13.4. \\ . However, observational evidence from the last decade has suggested 
that most MCs are l ikely to be at least moderately magnetical l y supe rcritical (e.g., 
Bourke et al. , 2001 ; Troland & Crutcher , 2008 ; Crutcher et al. , 2010l) . a conclusion 
also reached by theoretical arguments (see 35. 4. 21 and references therein). 

Another alternative was the proposal that MCs could be supported by turbulence, 
either hydrodynamical or MHD. However, since turbulence needs to be continu- 
ously driven, two variants have been considered for the driving: either it might be 
due to feedback from stellar sources internal to the clouds (cf. 35.4. Il l, or else to ex- 
ternal driving sources such as supernova shocks. However, as discussed in § 34.2.11 
and 15.4.21 the role of external turbulence seems more likely to be to drive MC for- 
mation, rather than to drive the strongly supersonic turbulence observed in GMCs, 
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because the turbulence induced in the forming clouds is only moderately supersonic, 
rather than strongly so (cf. ^6. It . 

The possibility of driving the turbulence by stellar feedback from i nside the 

cloud s has been extensively studied, both analytically and numerically (e.g.jNorman & Silk , 



1980 t jMcKeel Il989l: iLi & Nakamurai |2004 iNakamura & Lil l2007t ICarroll et al 



2009 , 2010t Wang et all 2010 ). In most such studies, it has been concluded that this 
feedback can maintain the clumps in near virial equilibrum. Studies of the SFR and 
the SFE under these conditions have often idealized the turbulence as being simply 
randomly driven, and have shown that in this case the SFE can be maintained at 
levels of a few percent, comparable to th e observed ones (e.g. jKlessen et al. . 2000; 
Vazquez-Semadeni et aUl2003a , l2005bl) . 



However, as discussed in the review by IVazquez-Semadeni (1201(1 numerical 
simulations of the momentum feedback from protostellar outflows have only con- 
sidered numerical boxes at the parsec (clump) scale, neglecting the infall from the 
clump's environment, which has been observed in GMC formation simulations. This 
adds a large amount of ram pressure to the system not included in those simulations. 
Thus, it seems that outflows cannot provide sufficient feedback to prevent the col- 
lapse of entire GMCs. 

The role of massive-star ioni zation feedback in the s uppo rt of GMCs has been in- 
vestigated semi-analytically bv lKrumholz et al. I (120061) and lGoldbaum et al.l (1201 lb 
considering the time-dependent virial theorem in the presence of feedback, and of 
feedback and infall, respectively, concluding that the clouds may oscillate around 
the virial equilibrium state for several Myr, until they are finally dispersed . How- 



ever, full numerical simulations of this problem ( Vazquez-Semadeni et al. , 2010h 



suggest that the infall is not suppressed, and instead that the regulation of the SFR 
occurs because most of the infalling material is evaporated before it can form furthe r 
stars, except in the case of the most massive (~ 10 6 M Q ) GMCs dDale et all 120121) . 



where supernova feedback may be also required to accomplish the dispersal of the 
clouds. 



Thus, it appears that the resolution of the IZuckerman & Palmer! d 19741) conun- 
drum lies not in the prevention of the global contraction of star-forming GMCs, but 
rather on the effect o f the feedba ck, and that this effect is, after all, essentially as 
initially suggested by lField dl970l) . The resolution of the objection by lMouschoviasi 
dl976l) . in turn, appears to lie in that the fra ction of mas s that is at very high den- 
sities (> 10 5 cm~ 3 ) is very small (see, e.g., BlitzL 1993 , sees. VII and IX), so that 
eventually HII regions may break out from the densest regions and ionize the rest of 
the MC. 



7 Summary and Conclusions 

In this contribution, we have briefly reviewed the role and interaction between the 
main physical processes present in the ISM: radiative heating and cooling, magnetic 
fields, self-gravity, and turbulence. The presence of radiative heating and cooling 
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implies in general that the gas behaves in a non-isentropic (i.e., non-adiabatic) way, 
and in particular it may become thermally unstable in certain regimes of density 
and temperature, where small (i.e., linear) perturbations can cause runaway heating 
or cooling of the gas that only stops when the gas exits that particular regime. This 
in turn causes the gas to avoid those unstable density and temperature ranges, and 
to settle in the stable ones, thus tending to segregate the gas into different phases 
of different densities and/or temperatures. In classical models of the ISM, only the 
stable phases were expected to exist in significant amounts. 

We then discussed some compressible MHD turbulence basics, and the produc- 
tion, nature and evolution of turbulent density fluctuations in polytropic (i.e., of the 
form P oc pTk) flows, discussing in particular the probability density function (PDF) 
of the density fluctuations, which takes a lognormal form in isothermal regimes, and 
develops power-law tails in polytropic ones. We also discussed the correlation (and, 
at low densities, lack thereof) between the magnetic field and the density as a con- 
sequence of the superposition of the different MHD wave modes, and the evolution 
of the mass-to-magnetic flux ratio (MFR) as density enhancements are assembled 
by turbulent fluctuations. 

We next discussed turbulence in the multi-phase ISM, noting that, since turbu- 
lence is an inherently mixing phenomenon, it opposes the segregating effect of ther- 
mal instability, causing the production of gas parcels in the classically forbidden 
unstable regimes, which may add up to nearly half the mass of the ISM, although 
the density PDF in general still exhibits some multimodality due to the gas' prefer- 
ence to settle in the stable regimes. The existence of gas in the unstable ranges has 
been established by various observational studies. 

Next, we discussed the nature of the turbulence in the different ranges of den- 
sity and temperature of the gas, noting that in the diffuse ionized regions, where the 
flow is transonic (i.e., with Mach numbers M s ~ 1), the gas appears to behave in an 
essentially incompressible way, exhibiting Kolmogorov scalings over many orders 
of magnitude in length scale. However, in the neutral atomic component, where the 
gas is thermally bistable, the flow is expected to exhibit large density and temper- 
ature fluctuations, by up to factors ~ 100, thus being highly fragmented. We also 
pointed out that large-scale compressions in the warm neutral gas, which may be 
triggered by either random turbulent motions, or by yet larger-scale instabilities, 
may nonlinearly induce the formation of large regions of dense, cold gas; much 
larger, in particular, than the most unstable scales of IT, which have sizes ~ 0. 1 pc, 
thus forming large cold atomic clouds that may be the precursors of giant molecular 
clouds (GMCs). This is because these clouds are expected to become molecular, 
gravitationally unstable, and magnetically supercritical at approximately the same 
time, so that when they reach a mostly molecular stage, they are likely to be under- 
going generalized gravitational contraction. 

The clouds are born internally turbulent and clumpy, and the resulting nonlin- 
ear density fluctuations ("clumps") may themselves become locally gravitationally 
unstable during the contraction of the whole large-scale cloud. Because they are 
denser, they have shorter free-fall times, and can complete their local collapses 
before the global one does, thus producing a regime of hierarchical gravitational 
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fragmentation, with small-scale, short-timescale collapses occurring within larger- 
scale, longer-timescale ones. It is thus quite likely that the flow regime in the dense 
molecular clouds corresponds to a dominant multi-scale gravitational contraction, 
with smaller-amplitude random (turbulent) motions superposed on it. 

The local collapses cause star formation (SF) that begins before the global col- 
lapse is concluded, and the ionizing feedback from the massive stars that form 
during this stage appears to be sufficient to erode and disperse the clouds before 
the entire mass of the cloud s is converted to stars, thus avoiding the objection by 



Zuckerman & Palmer! (1 19741) to free-falling GMCs, that they would form stars much 
more rapidly than observed. They do so, but only for short periods of time, before 
they get dispersed. 

We conclude that turbulence in the magnetized, multi-phase, self-gravitating ISM 
is an extremely rich and complex phenomenon, but whose (thermo)dynamics is be- 
ginning to be understood, together with its relation to the star formation process. 
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